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Abstract. We present a trajectory based interpretation for Young's experiment, 
the Arago-Fresnel laws and the Poisson-Arago spot. This approach is based on the 
equation of the trajectory associated with the quantum probability current density 
in the case of massive particles, and the Poynting vector for the electromagnetic 
field in the case of photons. Both the form and properties of the evaluated 
photon trajectories are in good agreement with the averaged trajectories of single 
photons observed recently in Young's experiment by Steinberg's group at the 
University of Toronto. In the case of the Arago-Fresnel laws for polarized light, the 
trajectory interpretation presented here differs from those interpretations based on 
the concept of "which- way" (or "which-slit" ) information and quantum erasure. 
More specifically, the observer's information about the slit that photons went 
through is not relevant to the existence of interference; what is relevant is the 
form of the electromagnetic energy density and its evolution, which will model 
consequently the distribution of trajectories and their topology. Finally, we also 
show that the distributions of end points of a large number of evaluated photon 
trajectories are in agreement with the distributions measured at the screen behind 
a circular disc, clearly giving rise to the Poisson-Arago spot. 



PACS numbers: 03.50.De, 03.65.Ta, 42. 25. Hz, 42.50.-p 
1. Introduction 

Very refined and ingenious interferometers and detectors for electrons [1], neutrons 
[2-4], atoms [5-7], molecules [7,8] and photons [9,10] have been devised to demonstrate 
that the quantum interference pattern can be build up by means of the accumulation 
of single detection events. Even before the realization of these experiments, De Broglie 
[11] and Bohm [12] argued that particles with mass possess simultaneously wave 
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and particle properties, and would move within an interferometer along trajectories 
determined by the guidance equation 

dr _ VS(r,t) 

dt m ' [ ' 

where S(r, t) is the phase of the particle wave function 

*(r,t) = \^(r,t)\e iS ^ t)/h (2) 

which satisfies the time-dependent Schrdinger equation. 

Using the method proposed by De Broglie and Bohm, Philipidis et al [13] plotted 
the trajectories of massive particles in the double slit experiment [13], Dewdney showed 
trajectories for neutrons inside a (neutron) interferometer [14], and Sanz and Miret- 
Arts explained the Talbot effect for atoms by plotting their associated trajectories 
behind a diffraction grating [15]. 

Starting from the works of Laukien [16], presented in [17], and Prosser [18,19], 
Davidovic et al explained [20] the emergence of interference patterns in experiments 
with photons by determining electromagnetic energy (EME) flow lines behind an 
interference grating. The equation of such EME flow lines reads as 

dr _S(r)_ 

ds cU(r) ' 1 ' 

where s denotes a certain arc- length along the corresponding path, S(r) is the real 
part of the complex-valued Poynting vector, 

S(r) = iRe[E(r)xH*(r)], (4) 

and U (r) is the time-averaged EME density, 

U(r) = \ [e E(r) • E*(r) + ^ H(r) • H*(r)] . (5) 

Here E(r) and H(r) are respectively the spatial part of the electric and magnetic field 
vectors, which satisfy Maxwell's equations and have been assumed to be harmonic, 
i.e., 



E(r) = E(r)e- iu;t 
H(r) = H(r)e 



(6) 



Davidovic et al pointed out [20] that it is useful to write the equation of the Bohmian 
trajectories for massive particles (pp) in terms of the probability current density, 

J(r, t) = — [tf W* - ^* V*] , (7) 

2im 

because from the latter form one can recast the guidance equation (pQ) as 
dr = J(r,t) 

dt |^(r,t)| 2 ' {) 

from which one may conclude that the equation for the EME flow lines and the 
equation of the Bohmian trajectories for massive particles have the same form. In 
other words, the Poynting vector in the case of photons plays the same role as the 
quantum probability current density in the case of particles with a mass. 

This analogy is even more apparent in cases where the spatial parts of the 
magnetic and electric fields can be expressed in terms of a scalar function that satisfies 
the Helmholtz equation [20,21]. 
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M. Gondran and A. Gondran [22] explained the appearance of the Poisson-Arago 
spot behind an illuminated circular disc using EME flow lines. In this way, they 
showed how such flow lines answer the question about diffraction phenomena presented 
two centuries ago by the French Academy " deduce by mathematical induction, the 
movements of the rays during their crossing near the bodies. " 

Recently, average trajectories of single photons in a double slit experiment 
were observed experimentally for the first time by Kocis et al [23]. Their result 
motivated us to apply the method of EME flow lines to numerically evaluate photon 
trajectories behind the double-slit grating with the same parameters as in the Kocis 
et al experiment. We show these results in Section O In Section [3] we study 
how polarizers put behind the slits affect the photon trajectories, thus providing a 
trajectory interpretation for the Arago-Fresnel laws. By adding orthogonal polarizers 
behind the slits, Kocis et al could observe average photon trajectories in the presence 
of orthogonal polarizers and check directly this interpretation. Section 2] is devoted to 
the trajectory based interpretation of the Poisson-Arago spot. 



2. EME flow lines - average photon trajectories in Young's interferometer 

Let us consider a monochromatic electromagnetic wave in vacuum incident onto a two 
slit grating located on the XY plane, at z = 0. In order to simplify the treatment we 
will assume that the electric and magnetic fields do not depend on the y coordinate. 
This assumption is justified when the slits are parallel to the y axis and their width 
along the y axis is much larger than the width along the x axis. In such a case from 
Maxwell's equations one obtains two independent sets of equations: one involving the 
H x and H z components of the magnetic field and the E y component of the electric 
field (commonly referred as ^-polarization), and another involving E x , E z and H y (H- 
polarization). As shown in [20], the electric and magnetic fields behind the grating 
are given by 

ifid^ i(3d^ 
E(r) = --^ + T ^e, + a*o y> (9) 

ia ia kf3e il ? 

uojiQ dz uj/^o dx ujiQ 
where ^ is a scalar function that satisfies the Helmholtz equation and the boundary 
conditions at the grating. The solution for \I/ can be written as a Fresnel-Kirchhoff 
integral, 

I T, poo 

9(X,Z) = J e ikz-in/4 ip(x',0 + )e ik(x - x ') 2/2z dx', (11) 

V ^7TZ J — oo 

where i(j(x\0 + ) is the wave function just behind the grating. 

We consider a grating with two Gaussian slits [23], so that the wave function just 
behind the grating is given by 

#r',0+) = Mx',0 + )+Mx',0 + ), (12) 



H ( r ) = — — e *"— — e * + ^— *<Hm ( 10 ) 



where 



V>iOr',0+)= [2^2) e-^'-^ 2 /^w(x' -fiua,), (13) 
M*, 0+) = (-^ V4 e-V-rt'^wtf - to, 02), (14) 
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Figure 1. EME density (a), phase of ^ (b) and photon trajectories (c) behind a 
two-slit Gaussian grating. The parameters are chosen as in the experiment carried 
out by Kocsis et al : a\ = a<2 = 0.3 mm, /ii = /i2 = 2.35 mm, a\ = a/2 = 1.8<ti 
and A = 943 nm. The initial polarization is lineal. The initial ^-coordinates 
for the trajectories are calculated from (|16|) . where u takes 19 equidistant values 
within the interval [0.02,0.98]. 



and w(x,a) is the window function, 

The EME flow lines (i.e., the average photon trajectories from the experiment carried 
out by Kocsis et al ) are obtained from equations (|3j)-((5j) and ([9])- ([T5]) . In figure [TJ 19 
photon trajectories per slit are shown. The initial x coordinates of the flow lines are 
chosen to be 

x s = Hi + (JiF~ x (u), (16) 

where i G [1,2] is the cardinal number of the slit and F _1 (i^) is the inverse of 
the Gaussian cumulative distribution function. If the variable u follows a uniform 
distribution, then the variable x s will have a Gaussian distribution with mean value 
\±i and variance 0{. 

By comparing the photon trajectories displayed in figure 1(c) with those 
experimentally inferred, shown in figure 3 in [23], we notice a very good agreement. 
The experimental average photon paths were reconstructed after performing a weak 
measurement on the momentum of an ensemble of photons and then a subsequent 
strong measurement of their position. In figure [2j we compare the experimental data 
coming from the measurement of the relative weak transverse momentum values, k x /k, 
as a function of the transverse coordinate at four different distances from the grating, 
with our theoretical curves obtained for three different window functions. 
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Figure 2. Transverse momentum along the transverse coordinate computed at 
four distances z from the two slits: (a) z = 3.2 m, (b) z = 4.5 m, (c) z = 5.6 m 
and (d) z = 7.7 m. The red solid line denotes the calculation with full Gaussians, 
while the blue and green lines refer to calculations where the outgoing beams 
were truncated Gaussians with a = 1.9a and a = 1.5a, respectively. To compare 
with, the experimental data (black circles) are also displayed. The parameters 
used for calculation are: a\ = 0.307 mm, ai — 0.301 mm, \i\ = 2.335 mm, 
fi2 = —2.355 mm, a\ = 1.5<xi, a2 = 1.5<72 and A = 943 nm. 



Since in our case the light propagates in vacuum, the Poynting vector can be 
identified with the density of electromagnetic momentum [24], so we have 

— = — . (17) 
k S K } 

The Fresnel-KirchhofT integral [17] can be integrated analytically for full Gaussians 

(for which the parameter of the window function a — >• oo), while the integration has 

to be done numerically for truncated Gaussians, as we did here. 



3. EME flow-line interpretation of the Arago-Fresnel laws 

According to a generalized version of the Arago-Fresnel laws, two beams with the same 
polarization state interfere with each other just as natural rays do, but no interference 
pattern will be observable if the two interfering beams are linearly polarized in 
orthogonal directions or elliptically polarized, with opposite handedness and mutually 
orthogonal major axes. The most direct way to verify these laws is by inserting 
mutually orthogonal polarizers behind the slits of a double-slit grating. 

The standard interpretation given to the disappearance of the interference fringes 
after inserting mutually orthogonal polarizers behind the slits is usually based on the 
Copenhagen notion of the external observer's knowledge (information) about the path 
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Figure 3. Photon trajectories behind a Gaussian double-slit grating followed by 
orthogonal polarizers. The parameters are the same as in figure ^ 



followed by the photon, i.e., the slit traversed by the photon in its way to the detection 
screen. Sanz et al [21] and Bozic et al [25] challenged this interpretation by explaining 
the first and second Arago-Fresnel laws considering EME flow lines behind the grating 
both in the presence and in the absence of polarizers. In both cases EME flow lines 
starting from slit 1 will end up in the side in front of slit 1, while those starting from 
slit 2 will end up in the side in front of slit 2, as also shown in quantum mechanics for 
matter particles [26]. However, the distribution of these EME flow lines is different 
in each case. In the absence of polarizers, the distribution shows interference fringes 
(figure 1(c) and figure 5 in [25]); in the presence of polarizers, the fringes are absent 
(figure 3 and figure 6 in [25]). 

As seen above, the average photon trajectories observed by Kocis et al. [23] in 
the absence of polarizers agree with our EME flow lines — the photon paths presented 
in figure [TJ In order to verify the interpretation of the Arago-Fresnel laws based on 
the EME flow lines, it would be interesting as well as challenging to experimentally 
determine the average photon paths in slightly modified experimental setup, by adding 
orthogonal polarizers behind the slits. In such a case, we expect that the corresponding 
experimentally inferred photon paths would look like the trajectories presented in 
figure [3l 

4. EME flow-line interpretation of the Poisson-Arago spot 

It is well known that the experimental observation of the so-called Poisson-Arago spotQ 
by Arago led to the acceptance of Fresnel's wave theory of light and the refutation 
of Newton's corpuscular theory of light. Now, by numerically evaluating EME flow 
lines behind a circular opaque disc, M. Gondran and A. Gondran found that these 
lines can reach the bright Poisson-Arago spot that appears at the center of the 
shadow region generated by such a disc [22]. These authors then argued that for 
a monochromatic wave in vacuum, the EME flow lines correspond to the diffracted 
rays of Newton's Opticks, thus concluding that after all Fresnel's wave theory may 
not be in contradiction with the corpuscular interpretation. This result also follows 

X This phenomenon is commonly regarded simply as Poisson spot. However, we have added the 
name of Arago in order to give him scientific credit, for it was him who provided the experimental 
evidence for this light phenomenon. 
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Figure 4. Photon trajectories in the XZ plane behind a circular disc in the 
XY plane, centered atx = y = z = and with radius R = 5 fim. The disc is 
illuminated by a monochromatic light with wavelength A = 500 nm. Because of 
the cylindrical symmetry of the problem, only trajectories having positive initial 
x coordinate are presented. 





Figure 5. Histogram of the end points of photon trajectories at the distance 
z = 15 mm behind a circular disc, having radius R = 0.5 mm, illuminated by the 
monochromatic light with wavelength A = 500 nm. 



from our evaluation of EME flow lines (figure HJ. Statistics of this lines (see figure [5]) 
agrees very well with the corresponding curve of light intensity behind the circular 
disc [22,27] determined by taking the square of the field function ^(P), which is given 
by the Rayleigh-Sommerfeld formula and Babinet's principle [22,28], 

* (*) = je ife * + jf ^ ( X " i) cos0dx Mdy M } , (18) 

where r = \/ (x — xm) 2 + (y — Vm) 2 + cos# = z/r, k = 2i\jr and integration is 
taken on the surface of the opaque disc S. 

Due to circular symmetry the integration over two variables in (fT8j) may be 
reduced to the integration over one variable. This simplifies and makes faster the 
numerical evaluation of the field function and the trajectories. In addition, this 
simplified formula makes possible the analysis of the dependence of the field function 
on the longitudinal z-coordinate. The details of this study constitute the subject of a 
forthcoming paper, where we will also present the Bohmian trajectories corresponding 
to a recent Poisson- Arago spot experiment performed with molecules [29] . 
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